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Recent experiments about the low temperature behaviour of a Single Wall Carbon Nanotube 
(SWCNT) showed typical Coulomb Blockade (CB) peaks in the zero bias conductance and allowed us 
to investigate the energy levels of interacting electrons. Other experiments confirmed the theoretical 
prediction about the crucial role which the long range nature of the Coulomb interaction plays in 
the correlated electronic transport through a SWCNT with two intramolecular tunneling barriers. 

In order to investigate the effects on low dimensional electron systems due to the range of electron 
electron repulsion, we introduce a model for the interaction which interpolates well between short 
and long range regimes. Our results could be compared with experimental data obtained in SWCNTs 
and with those obtained for an ideal vertical Quantum Dot (QD). 

For a better understanding of some experimental results we also discuss how defects and doping 
can break some symmetries of the bandstructure of a SWCNT. 

PACS numbers: 



I. INTRODUCTION 

In the last 20 years progresses in technology allowed 
the construction of several new devices in the range of 
nanometric dimensions. The known Moore prediction 
states that the silicon-data density on a chip doubles ev- 
ery 18 months. So we are going toward a new age when 
the devices in a computer will live in nanometer scale 
and will be ruled by the Quantum Mechanics laws. 

QDgiSi, which could play a central role within the 
quantum computation as quantum bits (QBIT's)^*^*^*^*! 
have been studied intensively in the last yearsS, thanks to 
the advances in semiconductor technologji^. 

QDs are small devices, usually formed in semiconduc- 
tor heterostructures, with perfectly defined shape and 
dimensiong^(< in diameter). They contain from one 
to a few thousand electrons and, because of the small vol- 
ume available, the electron energies are quantized. The 
QDs are useful to study a wide range of physical phe- 
nomena: from atomic like behaviouriiS*iiiiS to quantum 
chaos, to the Quantum Hall Effect (QHE) when a strong 
transverse magnetic field acts on the devic9i2*i^4i^*i&. 

Recently several scientists proposed a new carbon 
based technology against the usual silicon one. In 
this sense the discovery of carbon nanotubes (CNs) in 
1991— opened a new field of research in the physics at 
nanoscalesi^. 

Nanotubes are very intriguing systems and many ex- 
periments in the last decade have shown some of their 
interesting propertiesiS. An ideal SWCNT is a hexago- 
nal network of carbon atoms (graphene sheet) that has 
been rolled up to make a cylinder. The unique electronic 
properties of CNs are due to their diameter and chiral an- 
gle (helicity) parametrized by a roll-up (wrapping) vector 
(n, m). This vector corresponds to the periodic bound- 
ary conditionsSSiSi and gives us the dispersion relations 
of the one-dimensional bands, which link wavevector to 



energy, straightforwardly from the dispersion relation in 
a graphene sheet 

/ . , 7rm , 
e,n,k = ±7^1- 4cos(— ) cos( — ) + 4cos2( — ) 

(1) 

where iVf, is the number of periods of the hexagonal lat- 
tice around the compact dimension (y) of the cylinder. 
If the SWCNT is not excessively doped all the excita- 
tions of angular momentum m ^ (corresponding to the 
transverse motion ky = m/R) cost a huge energy of or- 
der (0.3 7 « 1 eV), so we may omit all transport bands 
except for the lowest one. From eq.(^ we obtain two 
linearly independent Fermi points ztK sa ^"^^ with a 
right- and a left-moving (r = R/L = ±) branch around 
each Fermi point (see Fig.(l.c)). These branches are 
highly linear in the Fermi velocity vp ~ 37/2« 8 x 10^ 
m/s up to energy scales E < D k, 1 eV. This linear dis- 
persion corresponds to that of a Luttinger Model for ID 
electron liquids. Many experiments demonstrated a LL 
behaviour'^^-^'^'^** in SWCNT^'- with a measurement of the 
linear temperature dependence of the resistance above a 
crossover temperature T^. 

In nanometric devices, when the thermal energy ksT 
is below the energy for adding an additional electron to 
the device (^jv = E{N) - E{N - 1)), low bias (small 
Vsd) transport is characterized by a current carried by 
successive discrete charging and discharging of the dot 
with just one electron. 

This phenomenon, known as single electron tunnel- 
ing (SET) or quantized charge transport, was observed 
in many experiments in vertical QDs at very small 
temperatureiSiiiiiS. In this regime the ground state en- 
ergy determines strongly the conductance and the pe- 
riod in Coulomb Oscillations (COs). COs correspond 
to the peaks observed in conductance as a function of 
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gate potential (V^) and are crudely described by the 
CB mechanismSi: the N — th conductance peak occurs 
wheni^ aeVg{N) ~ where a = ^ is the ratio of the 
gate capacitance to the total capacitance of the device. 
The peaks and their shape strongly depend on the tem- 
perature as explained by the Beenakker formula for the 
resonant tunneling conductancai^ 



"[ksT sinh(- 
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here /ii, ...,/iAr represent the positions of the peaks. 

Many experiments showed the peaks in the conduc- 
tance of QDs and a famous onc^^ also showed atomic-like 
properties of a vertical QD. There the "addition energy", 
needed to place an extra electron in a semiconductor QD, 
was defined analogously to the electron affinity for a real 
atom and was extracted from measurements. The typi- 
cal addition energy ranges from 10 to 20meV, while the 
disappearance of the COs happens above « 50 — K. 

The transport in CNs often differs from the quantum 
CB theory for QDs, because of the one-dimensional na- 
ture of the correlated electrons: so we could need a pe- 
culiar theory for resonant tunneling in LLs^^^^% More 
recently a novel tunneling mechanism^S was introduced 
i.e. correlated sequential tunneling (CST), originating 
from the finite range nature of the Coulomb interaction 
in SWNTs, in order to replace conventional uncorrelated 
sequential tunneling. It dominates resonant transport 
at low temperatures and strong interactions and its pre- 
diction agree with experiments^!. In the experiment a 
short nanotube segment was created with an addition 
energy larger than the thermal energy at room tempera- 
ture (Tfl), so that the SET can be observed also at Tr. 
The conductance was observed to follow a clear power 
law dependence with decreasing temperature, pointing 
at a LL behaviour in agreement with Ref."^*^. 

An interesting observation about the transport in CNs 
concerns the long-range nature of the Coulomb inter- 
action, which induces dipole-dipole correlations between 
the tunneling events across the left and right barrier'^''. 
The crucial question of the range of the interaction in 
CNs was investigated in different transport regimes e.g. 
in order to explain the LL behaviour of large Multi Wall^S 
and doped^^ CNs. 

However many experiments showed COs: e. g. in 1997 
Bockrath and coworkers^'* in a rope of CNs below about 
10 ~ K observed dramatic peaks in the conductance as 
a function of the gate voltage according the theory of 
single-electron charging and resonant tunneling through 
the quantized energy levels of the nanotubes composing 
the rope^. In this regime also a SWCNT behaves as an 
artificial atom and reveals its shell structur0^(the data 
were taken at 5mK). Recent measurements report clean 
closed nanotube dots showing complete CB^S, which en- 
able us to deduce some properties from the addition en- 



ergy of SWCNT and discuss the role which the Coulomb 
interaction could play in a ID system at small tempera- 
tures (T = 0.1-^0.3-/0- 

In this paper we analyze the effects of a long range 
electron electron interaction in order to determine the 
addition energy for models chosen by a vertical QD and 
a SWCNT in the Hartree Fock (HF) approximation. 

In section II we introduce microscopic models and cor- 
responding Hamiltonians for SWCNTs and QDs also fo- 
cusing on a theoretical model for the interaction potential 
which interpolates between short and long range type in- 
teraction. In section III we show our results about the 
SWCNT and QDs and discuss the effects of asymmetries 
experimentally observed^^. 



II. MICROSCOPIC APPROACH 

As showed in eq.(^ near each Fermi point we ob- 
tain that the CN could be represented by a typical Lut- 
tinger Hamiltonian with linear branches depending on 
k — k + aKp {a — ±1 labels the Fermi point). We in- 
troduce the operators that create the electrons near one 
of the Fermi points (label a) belonging to one of the two 
branches (label C corresponding to the sign of A;) , ^ 

and - . In terms of these operators the free and in- 

a,k,s 

teraction Hamiltonians can be written as 



(3) 



Hint — (^k,p],s'i'l)'^lii,k+q,s'^2,p-q,s''^a3,P,s''^a4,,k,i^) 
{ai} k.k' ,q,s,s' 

The interaction V^^^,'^(g) plays a central role in order 
to determine the properties of the electron liquid. 

A very crucial question is the effective range of the 
potential and its possible screening in a CN. If we denote 
by X the longitudinal direction of the tube and y the 
wrapped one, the single particle wave function for each 
electron reads 



V2^ 



where Ui^^a{x,y) is the appropriate linear combination of 
the sublattice states p = ± and L the length of the tube. 
So we can obtain a simple ID interaction potential as 
follows: 



'lf^}{x-x')= I dydy' ul^^{x,y)ul,^^{x',y') 



xU^-'Ax - x' ,y~y')u^,^^{x ,y')u(^^^{x,y) . 

These potentials only depend 0T^ x — x' and the ID 
fermion quantum numbers while C/q '' ~ x',y — y') is 
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obtained from a linear combination of U{x — x' ,y — y' + 
pdSp _pi) sublattice interactions'^^. 

Because of the screening of the interaction in CNs and 
the divergence due to the long range Coulomb interac- 
tion in ID electron systems, it is customary to introduce 
models, in order to describe the electron electron repul- 
sion. The usual model is the so called Luttinger model, 
where the electron electron repulsion is assumed to be 
a constant in the space of momenta, corresponding to a 
very short range ID potential (Dirac delta). In order to 
analyze the effects of long or short range interactions, we 
introduce a model for the electron electron potential de- 
pending on a parameter r, which measures the range of 
a non singular interaction; it has as limits the very short 
range potential (r — > 0, delta function) and the infinite 
long range one (r — > oo, constant interaction). So we can 
conclude that our general interaction model ranges from 
the very short range one to the infinity long range one 
and eliminates the divergence of the Coulomb repulsion. 
In this sense we suppose that our model is good for de- 
scribing the interaction, if we do not take in account the 
IR and the UV divergences. 



Ur{\x^x'\)=Uo 



2r 



(5) 



The interaction between two different electrons with mo- 
menta k and q follows from the previous formula 



V{p^\k~q\ 

In the limit L — )■ c 
form of eq. ^ 



1 

l2 



dx 



e'PyU{y)dy 



we can calculate the Fourier trans- 



(6) 



The scattering processes are usually classified accord- 
ing to the different electrons involved and the coupling 
strengths g are often taken as constants. This assump- 
tion corresponds to the usual Luttinger model, so we fol- 
low this historical scheme, in order to classify the in- 
teractions. The backscattering g^''^ involves electrons 
in opposite branches with a large momentum transfer 
{q « 2kp) so gl''^ « V{2kF). The forward scattering 
occurs between electrons in opposite branches 172 with a 
small momentum transfer {q ^ 2kp) so 52'* ~ ^(0) (or 
g^'^ « V{2tt/L)). The forward scattering in the same 
branch (54) involves the pairs (fc ft! kpjP kp) and gives 
gl » V{0) - V{p - k) and .gf » V{Q). Further below, 
where we discuss the electron electron interaction, we re- 
call the values of the g constants. 

Now we want to introduce the analogous model for 
a semiconductor ideal QD. Usually we describe the dot 
like a 2D system with an harmonic confinement poten- 
tial V{r) = jTTi* uJd^r^ according to measurements^^ that 



demonstrated that vertical QDs have the shape of a disk 
where the lateral confining potential has a cylindrical 
symmetry with a rather soft boundary profile. Under this 
hypothesis the quantum single particle levels depend just 
on n = n+-|-n_ (the angular momentum is m = n+ — n_) 

En = hcjd{n + i). 

The symmetry leads to sets of degenerate single-particle 
states which form a shell structure: each shell (£„) has 
2{n + 1) degenerate states so that the shells are com- 
pletely filled for = 2, 6, 12, 20, etc. electrons in the dot 
(Magic Numbers). 

The many body Hamiltonian corresponding to eq.Q 
and eq.Q) has the form 



H 



5 E 

Q,/3,7,(5 



(7) 



Here a = (n, m, s) denotes the single particle state in 
the single particle energy level Sa, c^. creates a particle 
in the state a and ha = cj^Cc is the occupation number 
operator. In the following sections we discuss the essen- 
tial question regarding the electron electron interaction 
in the dot {V^'^,) and analyze in detail the screening of 
the effective potential. 



III. SINGLE WALL CARBON NANOTUBE: 
LOW TEMPERATURE BEHAVIOUR AND 
COULOMB BLOCKADE 

Before proceeding with the calculations, we want to 
point out that the real band structures of measured CNs 
show some differences with respect to the ideal case dis- 
cussed (eq.©): in order to clarify this point we shortly 
discuss the model and the results of two recent experi- 
ments. 

To begin with, we have to introduce a quantization 
due to the finite longitudinal size of the tube (L) in the 
dispersion relation eq. . The longitudinal quantization 
introduces a parameter which also gives a thermal limit 
for the Atomic Like behaviour: in fact k wavevectors 
have to be taken as a continuum if KgT is as a critical 
value Ec = vpj^ and as a discrete set if the temperature 
is below (or near) Ec- 

After the quantization we obtain shells with an 8-fold 
degeneracy due to a (spin symmetry), a {K,—K lattice 
symmetry), C ((fc — K), {K — k) Luttinger symmetry). 

Recent experiments do not support such a high sym- 
metry and different hypotheses were formulated in order 
to explain this discrepancy. 

According to Cobden and Nygard2& "the sole orbital 
symmetry is a two-fold one, corresponding to a K-K' sub- 
band degeneracy and resulting from the equivalence of the 
two atoms in the primitive cell of graphene structure". 
Experimentally one can answer this question by observ- 
ing the grouping of the peaks in plots of the conductance 
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FIG. 1: The dispersion relation and the quantized levels. The 
boxes in figure represent energy levels and can be filled by a 
pair of electrons with opposite spins, a) The general case 
without any degeneracy, b) The 4-fold degeneracy case with 
SsM ~ 0. c) The 8-fold degeneracy case, d) Differences in the 
splitting due to the comparison between Ssm and Ae. 



versus the gate potential. However in the experiment no 
four-fold grouping was observed because degeneracy was 
lifted by a mixing between states due either to defects or 
to the contacts. 

A different experiment^ displays conductance peaks 
in clusters of four, indicating that there is a four fold de- 
generacy. In refjSS: two different shell filling models are 
put forward: the first one, when the subband mismatch 
dominates, predicts that the spin in the SWCNT oscil- 
lates between S = and S" = 1/2. 

In order to take into account the strong asymmetries 
measured experimentally we modify the dispersion rela- 
tion. A first correction has to be introduced because of 
the "longitudinal incommensurability": in general K is 
not a multiple of n/ L so K = {N + SN)^ with 6N < 1 
and the energy shift is Ae — vp^^^. A second correc- 
tion is due to the subband mismatch (Ssm)- The single 
electron energy levels are 



£La,p = hvplly +pK\ 



JSM 



(8) 



where p = ±1. 

Each choice of parameters gives a different degener- 
acy for the quantum levels: the 8-fold degeneracy corre- 
sponds to Ssm = and K = nj^; the 4- fold degeneracy 
is found if we put just Ssm — and the 2 fold degeneracy 
represents the general case (see Fig. (I)). 



A. High band structure symmetry: damping in the 
addition energy oscillations 



laws of our Hamiltonian and find the Number of electrons 
(N), the Energy, the total linear momentum K ~ and 
the spin S, Sz — 0- 

The Fermi sea corresponds to the state where all 
the shells with energy below the Fermi energy [Ep = 
VFnph/L) are totally filled {Np = 4np). This state will 
be our ground state 'i>o{Np) and this is true also for in- 
teracting electrons in absence of correlation (i.e. in the 
HF approximation ). 

The effects of correlation will be discussed in a further 
article, here we have to explain the range of validity of 
our approximation. The first thing to consider is the 
interaction strength g « ^^(0) — V{2kp) compared to the 
kinetic energy vph/L. The HF approximation is valid 
if g <C Vph/L. However also the temperature plays a 
central role, in fact if the temperature increases we have 
to take into account more excited states (a sort of thermal 
cut-off corresponds to the energy kgT) so if we are at a 
very low temperature we can assume the Fermi sea state 
as the ground state. 

At this point we are able to calculate the addition 
energy (E-^ ) following the Aufbau sequence explained 
in Appendix A. E-^ has a maximum for some num- 
bers 4,8, 12, .An (n integer) due to the shell filling. The 
shells are filled sequentially and Hund's rule determines 
whether a spin-down or a spin-up electron is added so 
that the singlet (5 = 0) energy for a, 4n + 2 system is 
always greater than the triplet {S = 1) one. Obviously 
this is an effect of interaction and is quite different for 
the long and short range models. 

As we show in Fig. (2) the oscillations due to Hund's 
rule correspond to the short range potential while the at- 
tenuation of these oscillations when the number of elec- 
trons in the ID system increases is due to the long range 
interaction. So we can draw the following conclusions: 

o The 4— fold degenerate model predicts oscillations in 
the addition energy due to the Hund's rule quite similar 
to the ones observed in QDs. 

o The oscillations periodicity is 4 for this model (8 for 
a system with two Fermi points) 

o The oscillations amplitude is due to an exchange term 
(proportional to the short range interaction). 

o The effect of a long range interaction is a damping 
of the oscillations when the number of electrons in the 
system increases. 

The model with 8-fold degeneracy (see Fig.(l)) has two 
basic symmetries: k —k, KL/tt — Nk and usually the 
interaction between the electrons with momenta near K 
and the ones with momenta near —K is very small so that 
we have two independent 4-fold degenerate Hamiltonians 
{np ^ N/4« Nk). 



In the discussion which follows we analyze the effects 
of the range of the Coulomb interaction in a simplified 
system with just two linear symmetric branches. In this 
model each shell is filled by 4 electrons with opposite 
momenta and spin. We can look for the conservation 



B. Asymmetric model 

Now we have to analyze models without symmetries 
by using eq.® in the HF approximation. 



FIG. 2: On the right we show analytical Aufbau results for 
the addition energy versus the number of electrons of a 4- 
fold degeneracy model corresponding to different values of 
the range r (r = dark gray dashed line, r = 0.4 black line, 
r — 1 gray line). We show how the damping in the oscillations 
is due to a long range interaction while it does not appear for 
a r = model. Our predictions can be compared with the 
measured addition energy in the Cobden Nygard experiment 
displayed on the left. 



In order to introduce the electron electron interaction, 
we take in account just two (f/" and g^) of the many 
constants that we introduced in section II from cq.jnj. 
For allowing to compare easily our results with those in 
ref.^^s as well as with experiments, we give our results in 
terms of Vq and J, obtained as a linear combination of 
the g constants. 

Following the usual method, in order to calculate the 
energy levels, we put 



^Vn-J 



A = 



vph 
~2L' 



As = SN 



Vph 



Here Ae is the incommensurability shift (Ae < 0.5A). 
The single particle energies have a different structure for 
SsM > 2A£ and 5sm < 2Ae, as we show in Fig.(l.d). 
From the experimental data^Si we obtain 

Vo~U + 6U + Jexp ~ -42 and J w Jexp - 2SU w .05 

where we assume U = .22, 6U — .05 and Jexp = .15 in 
units of A. Under these conditions J is always less than 
the level spacing. 

Now we can write the Hamiltonian of the nanotube 
depending on these parameters'^ 



(9) 



+ Vo—^ --J 2^ > 



Since J is less than the level spacing the energy that we 
need, in order to add one electron is e, corresponding to 
the lowest empty energy level with an interaction energy 




FIG. 3: The asymmetric model calculation for the COs (con- 
ductance vs. gate voltage) at different temperatures expressed 
in terms of A calculated following the classical CB theory. 
Theoretical calculations show that the fine structure peaks 
are appreciable just for very low temperatures. 



Starting from the Hamiltonian eq.@ we are able to cal- 
culate the ground states of the many electron system for 
various N . 

Our results can be compared to the experimental re- 
sults where strong asymmetries were found. In Fig. (3) we 
show the peaks corresponding to an asymmetric model 
to which we apply the classical theory of CB, i.e. eq.ljSJ). 
The fine structure with 8 periodicity is destroyed by ther- 
mal effects: it is appreciable at smaller temperatures (we 
assume about T, w 300-^500mi4r for a Nanotube's length 
of about 100 SOOnm) and disappears at a temperature 
T « 4Ts (T « 1.2 ^ 2.QK see ref,^) where just a 4 pe- 
riodicity appears. So, we conclude that we should not 
be able to observe any small asymmetries effects if the 
temperature increases. Fig. (3) shows how the Coulomb 
peaks in the conductance disappear when the tempera- 
ture increases. In a future article we will show how the 
end of the CB regime corresponds to the beginning of 
another one. 

Now we want point out the limits of the approach used 
in the previous sections: the HF approximation ignores 
the effects of correlation. This corresponds to the Fermi 
Liquid theory and gives good results just if we can assume 
that the interaction is smaller than the kinetic energy 
{g < Tivp)- 

Thermal effects are quite important too. In fact the 
temperature appears in the Beennaker formula and is 
responsible of the disappearance of the Coulomb peaks. 
However, if the temperature is higher for a SWCNT, it is 
the same Beennaker formula which fails, because the HF 
calculated energy levels are very different from the real 
energy levels of the electron system. When T is above a 
critical value 



T, = 



Lkn 



(10) 



we cannot consider the Fermi sea as the ground state 
of the electron system in a SWCNT because some other 
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states with the same linear momentum K = and differ- 
ent kinetic energy are also available for the system. So 
for temperatures above the critical value we have to take 
into account strong effects of correlation. 



C. Quantum Dots and Long Range Interactions 

As we did for a SWCNT we calculate the addition en- 
ergy of a QD in the HF approximation. 

As we discussed in section II peaks and oscillations in 
the measured addition energy correspond to the ones of a 
shell structure for a two-dimensional harmonic potential. 
However, in experiments^^, high values of the addition 
energy are observed also for TV = 4, 9, 16, etc. correspond- 
ing to those values of the Number of electrons in the 
dot for which, respectively, the second, third and fourth 
shells are half filled with parallel spins in accordance with 
Hund's rule. Half filled shells correspond to a maximum 
spin state, which has a relatively low energjiifl*ii. 

We compare a non-interacting model with models that 
include Coulomb interactions especially the exchange 
term. In Fig. (4) we plot the calculated addition energy 
as a function of the number of electrons in the following 
three different cases: 

o The non-interacting model gives us just the peaks 
corresponding to the Magic Numbers. 

o The model with constant parameters (5" and g^) 
corresponds to a short range interaction (Dirac S). The 
Coulomb exchange term in HF gives the Hund's rule and 
allows us to explain the oscillations in the addition energy 
as we show in Fig. (4. a). 

o The third model is the long range interaction one, 
where we introduce two measured effects that are both 
due to the long range Coulomb interaction. The first one 
is due to the classical capacitive effect shown in Fig.(4.b 
and c) as a continuous line. The second one is due to the 
long range correction to Coulomb exchange. 

We could compare these results to the experiments 
(e.g. see reSii). 



Model: constant interaction 
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FIG. 4: In a),b),c) results from three different models of in- 
teracting electrons in QDs are displayed. The simple HF cal- 
culation in a) has to be corrected because it does not take 
into account the classical capacitive effect. In order to deter- 
mine it we recall that the electrons in the dot give a charge 
droplet of radius Kd not fixed, so that we can approximate it 
with a disk of capacitance C = CoRd- The value of Rd can 
be calculated from classical equations, Rd oc + 1. So we 
add the classical term to the damped oscillations due to long 
range affected Coulomb exchange and obtain the b) and c) 
plot in the figure. 



So we have to limit ourselves to describing an electron 
system only when it is uncorrelated. This is true only 
if the interaction strength is small {g << vph/L)'^^ and 
the temperature is below the critical value Tc (see ea[TU)l . 
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IV. CONCLUSIONS 

In this work we have analyzed some properties of a 
ID (SWCNT) and 0D-2D (QDs) electron systems by in- 
troducing a model of interaction capable of interpolating 
from short to long range, in order to analyze the effects 
of the long range component of electron electron interac- 
tions. 

We have obtained that a damping in the oscillations of 
the addition energy could be predicted for these models 
corresponding to the presence of a long range interaction. 
The results can be compared to the experiments and we 
discussed the limits for the experimental observation of 
our predictions in SWCNT. In fact the breakdown of the 
Fermi Liquid opens different regimes where the tunneling 
has to be considered as resonant tunneling in LLs22j^. 



APPENDIX A: AUFBAU AND ENERGIES FOR 
THE 4-FOLD DEGENERATE MODEL 

We start by taking into account that each electron in 
the |fc| > kp state interacts with all the electrons in the 
filled shells with k below kp so that we can have the 
following two (or four) "Ehf terms: 

p—i P——1 

If we introduce np = —kp and consider the direct term 
of the interaction, we conclude that 
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2nF- ^ J\p-k\{^ - ^pfi) 

Obviously if wc add an electron in the lowest 
energy empty shell we obtain Y,{kF,kF + x) ~ 

Vo(AnF-Y.liT'J.)- 

For each shell we can consider the internal interaction 
energy as 

W{k) = 2^0(3 - J2k) W{kF) = J2 ^(P) 



where /Ujv+i = -Bjv+i — and 7 : 



J2kp 
Uo 



The addition energy is also calculated as 
Mjv+2 — Mjv+1 so that 



so that the total energy of a system with uf filled shells 
reads 

Eir,p=2VFnF{nF + l)+W{kF)+^ X] 

a. Energy levels 
When we add one by one the electrons to the Auf- 



(^kF + j)+^{kF,kF + j) 
[kp + ^)+ ^{kp, kF + 2) + Uo{l - 7) 
[kF + 2)+ '^F + J^)+ 2C^o 

(kF + J) + ^{kp, kF + 2) + f/o(3 - 7) 
(kp + 2^) + T,{kF + ^,kF + 2^) 



system we obtain 


fJ-inp+l 


= hvp I 




= hvp i 




= hvp ( 




= %Vf ( 


M4nj.+5 


= hvp ( 



<P + 1 = C^0(1 - 7) <,+2 = t^0(l + 7) 



i^4l.+3 = C^o(l-7) 



C/o(l + 7) + 



hvpTT 

L 
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